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Abstract 

By a suitable shifting-the-mean parametrization at the Dirichlet series level and Delange's Taube- 
rian theorems, we show that the number of factors in random ordered factorizations of integers is 
■ asymptotically normally distributed. 

t-H ' 

^ ■ 1 Introduction 

(N ■ 
O ' 

Q\ \ Let 8? be a fixed subset of {2,3, . . . }. Let a(n) denote the number of different ways of writing n 
O ■ as the product of ordered sequences (ni, . . . ,n k ) of integers in Define a(l) = 1. Let A(N) : = 
^ • ^2 1<n<N a(n). Assume that all A(N) factorizations of an integer < iV are equally likely; denote by Yn 
^ . the random variable counting the number of factors in a random factorization. We prove in this paper that 
the distribution of Y N is asymptotically normal under very general conditions on 
Denote by Vis) the Dirichlet series of & 



V{s) := £ 



n 



Assume the abscissa of convergence of V(s) is k. Then k = — oo if | &\ < oo and < k < 1 if | = oo. 
Note that V{k) < oo. 

Our main result is as follows. 



*Part of this work was done while both authors were visiting Institut Mittag-Leffler, Djursholm, Sweden. They thank the 
Institute for hospitality and support. 
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Theorem 1. Assume \ > 2 and 1 < V(k) < oo; thus there exists p > max{K, 0} such that V(p) = 1. 
Let ^ 

v := ~v\p)' and ° 2 := ^ V "^ ~ M " 



Then 



where — > stands for convergence in distribution and tsV(0, 1) denotes the standard normal distribution. 
The mean and variance of Yn are asymptotic to \i log A" and a 2 log N, respectively. 

By Cauchy-Schwarz inequality, a 2 > 0. We will indeed prove convergence of all moments. 

The case when \&>\ = 1, say & = {d}, d > 2, is exceptional. In this case, a(n) = 1 when n is a 
power of d, and a(n) = otherwise. Then is uniformly distributed on the integers {0, . . . , [\og d N\}, 
and thus Y N / log iV converges in distribution to a uniform distribution on [0, 1/ log d] , and therefore, Y N is 
not asymptotically normal; further, the moment asymptotics differ from those in the theorem. 

Ordered factorization problems in connection with that studied in this paper have a long history, tracing 
back to at least MacMahon's work (see lfT9lO in the early 1890's; later they were in most publications 
referred to as Kalmar's problem of "factorisatio numerorum" (see lfT51 [T3l0 . Diverse properties of such 
factorizations have then been widely investigated, often in quite different contexts, one reason being that 
ordered factorizations are naturally encountered in many enumeration problems. For example, when & = 
{2, 3, . . . }, A(N) + 1 equals the permanent of the Redheffer matrix; see [|23l ; also they appeared as the 
lower bound of certain biological problems; see [|2TT|. See also [E1[T7]| for more information and references. 

The first paper dealing with general ordered factorizations beyond the subset = {2, 3, . . . } similar 
to our setting was Erdos [6J, extending previous results by Hille [12); see also lfT6l . Asymptotic normality 
of the special case of Theorem Q] when = {2, 3, . . . } was treated in lfl4ll . In this case, p ta 1.7286 being 
the unique root > 1 that solves ((s) = 2, where £ denotes Riemann's zeta function. The proof given there 
relies on the determination of a zero-free region of the function 1 — z(((s) — 1), which in turn involves 
deep estimates from trigonometric sums (see also lfT8lO . Such refined estimates are for general & hard to 
establish. We replace this estimate by applying purely Tauberian arguments of Delange (see [4]), which 
require only analytic information of the involved Dirichlet series on its half-plane of convergence. 

We will use Dirichlet series and the method of moments and derive asymptotic estimates for central 
moments of integral orders, which, by Frechet-Shohat's moment convergence theorem, will suffice to 
prove the theorem. 

Proposition 1. For k > 

k ( fc 1 

lim E ( Yn -t^l N ) = niSeVen CD 

v ^v^at j y y if k is odd, 

A straightforward application of the Tauberian theorem does not provide precise asymptotics for central 
moments beyond the first due to cancellation of major dominant terms and due to the fact that no error term 
is generally available via application of Tauberian arguments. The new idea we introduce in this paper is to 
take into account the feature that the mean is logarithmic and to shift-the-mean on the associated Dirichlet 
series, which nicely incorporates the cancellations of higher central moments in a surprisingly neat way. 
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It thus avoids completely the messy calculations and cancellations that the usual method of moments 
faces when dealing with higher central moments. A similar idea was previously applied to characterize 
the phase change of random m-ary search trees, where a nonlinear differential equation with an additive 
nature was encountered; see Q, However, the tools used there are complex- analytic, in contrast to the 
purely Tauberian ones used here. 

In the special case when & = {2,3,...}, our result may be interpreted as saying that the property of 
the zero-free region for the Dirichlet series 1 — z(((s) — 1) lies much deeper than the asymptotic normality 
of the random variables Y^. 

In addition to the number-theoretic interest per se of our results, we believe that the approach we use 
here also offers important methodological value for the study of similar problems. In particular, not only is 
the use of the method of moments very simple, but no analytic properties of the Dirichlet series beyond the 
abscissa of convergence are needed, which largely simplifies the analysis in many situations. For example, 
our approach can be readily applied to other factorizations such as branching or cyclic factorizations with 
algebraic or logarithmic singularities (see lfT3lO . It can also be extended, coupling with suitable Tauberian 
theorems, to deal with ordered factorizations in additive arithmetic semigroups (see lfT6lO and components 
counts in ordered combinatorial structures (see [8]). Another possible extension is to the analysis of 
Euclidean algorithms; see IfTI fTOllZZll . 

In the periodic case when & C {d k : k > 1} for some d > 2, V(s) has period 2iri/ logrf and the 
usual Tauberian theorem does not apply. Instead we write V(s) = V(d~ s ), where V(z) := J2d k e,$> zh > an ^ 
transform the multiplicative nature of the problem into an additive one on random compositions by taking 
logarithms. We replace the Tauberian theorem by the singularity analysis of Flajolet and Odlyzko (see 
or [|9l Chapter VI]) in the proof; the details are similar to the proof below, but simpler, so we omit them. 
In the rest of the paper we thus assume, for every d > 2, 

9 % {d%>i. (2) 

2 Dirichlet series, Delange's Tauberian theorem, and proofs 

Generating functions. Let a m (n) denote the number of ordered factorizations of n into exactly m fac- 
tors. Then 

P(s) m = J2a m (n)n- s , 

n>l 

in formal power series sense; analytically, we can take s satisfying V($t(s)) < oo. Thus if e lt ^V(^:(s)) < 
1, then by absolute convergence 

J2n- s J2a m (n)e- = £ e™^)™ = j _ ' • (3) 

n>l m>0 m>0 ^ ' 

Delange's Tauberian theorem. We need the following form of Delange's Tauberian theorem (see flU or 
iH Ch. Ill, Sec. 3]). 

Let F(s) := J2 n >i d(n)n~ s be a Dirichlet series with nonnegative coefficients and convergent 
for 3? (s) > q > . Assume (i) F(s) is analytic for all points on 3?(s) = q except at s = q; (ii) 
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F(s) = + H(s) 09 >0), 



for s ~ 9ft(s) > g, 

G{s 
[s-q) 

where G and H are analytic at s = g with G(g) ^ 0. Then 

G(g) 



E^-^^OogiV)'- 1 . (4) 



n<JV 



Asymptotics of A(iV). Taking 2 = in ©, we obtain the Dirichlet series for a(n) = Xlm>o a ™( n ) 

n>l \ ' 

as long as 3?(s) > p. Note that the non-periodicity assumption © implies that 'P(s) 7^ 1 for all s with 
3?(s) = p but s ^ p. Hence ^4(s) is not only analytic in the open half-plane {s : 3?(s) > p} but also on 
the boundary {s : 3?(s) = p} except at s = p. The same holds true for all Dirichlet series we consider 
below. 

Now for our V(s), since V'(p) = — J2 n &^ n ~ P ^°S n < 0> we see tnat ^P{ s ) nas a simple zero at s — p, 
and thus A(s) has a simple pole at p with 

n _ 1 -1 
as s — > p. Hence Delange's Tauberian theorem applies and we obtain 

Furthermore, we also have 

^a(n) (log n) fc ~ P N p {\ogN) k ~ A{N){\ogN)\ 

n<N P 

either by repeating the same procedure for the Dirichlet series 

= J>(n)(logn) fc n-* = , (6) 

n>l ^ ' 

or by using directly ©. The estimate will be used later. 

The expected value of Y^. By taking the derivative with respect to z on both sides of ©> we obtain 

V{s) 



ma m (n 



n>l m>0 v v " 

Delange's Tauberian conditions being easily checked as above, we then obtain 

E ^ = AW) E J2 mam ^ ~ ^ogN. 

^ ' n<N m>0 
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Shifting-the-mean at the Dirichlet series level. For higher central moments, the idea we will use can 
formally be described by using Perron's integral representation as follows (using © and for simplicity 
assuming temporarily that N is not an integer). 

i i>c+ioo jiTs—uz I 

ds, 



2itiA{N) J c _ ioo s (1 +fiz/s)(l - e z V(s + fiz)) 

where c is suitably chosen; the fact that the mean being of order log N is crucial here. We then formally 
expect that 



U rc+ioo yys 



where Qfc(s) is the coefficient of z k in the Taylor expansion (in z) of 

1 



(l + fiz/s)(l -e z V(s + fj,z)) 

While all steps can be easily justified (as done below), we cannot directly apply Delange's Tauberian 
theorem to Qk{s) here since each Qk (except Qo(s)) is not a proper Dirichlet series, but involves additional 
powers of This can be resolved as follows. 

Shifting-the-mean at the coefficients level. We look at the "translation" of the preceding parameter- 
shift at the coefficient level. By definition 

A(N)E(e (YN -^ N)z ) = J2J2 am ^ e{m ~^° gN)Z 

n<N m>0 

_ am ^ n ^ e (m-^logn)z-iiz\og(N/n) 
n<Nm>0 

Let 

h(n) := ^ a m{n){m - [L log nf . 

m>0 

Then, by taking the coefficients of z k on both sides, we obtain 

A(N)E(Y N -^ogN) k = 0)(-/^E^)(log^V ■ (8) 

0<i<k ^ ' n<N ^ ' 

We will see that the growth order of J2 n<N bk(n) is the power N p times an additional logarithmic term; 
it then follows that the weighted sum on the right-hand side is of the same order by a simple partial 
summation (see below for more details). 

Now observe that (assuming again that N is not an integer) 

where c is taken to be any real number greater than the abscissa of absolute convergence of the function 
defined by the series ^ 7 >i a U)j s - So, this, together with ©, justifies ©. 
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A probabilistic interpretation. Given N, consider a random factorization of a number n < N, namely, 
a random product p\ ■ • -p m < n with all pj e (uniformly distributed over all A(N) possible factoriza- 
tions). Let Yjv be the number of factors (= m) and un be their product (= n). Then 

A(N)E(Y N - fi\ogu N ) k = ^6 fc (n), 
which gives a probabilistic interpretation of the partial sum. 



The Dirichlet series of b k {n). Define the Dirichlet series 

M k (s) := ^2b k {n)n' s = ^n~ s a„(m)(m - /i log n) fc . (9) 

n>l n>l m>0 

Note that a m (n) > implies that m < log 2 n, so that 

(m - n\ogn) k = O ((logn) fc ) , 

for all non-zero terms. Hence Ai k (s) is absolutely convergent when 3?(s) > p because A(s) = 1/(1 — 
V(s)) is. Now if 3?(s) > p and |z| is sufficiently small, then, by ©, 

E = E „-, Eo „ (m)e <™-^ = 1 (10 ) 

fe>0 n>l m v ^ y 

With these M k (s), the generating function Q k (s) can be decomposed as 

0<£<fc v 7 

Our strategy will then to apply Delange's Tauberian theorem to Ai k for even and some auxiliary 
Dirichlet series for odd k, and then the asymptotics of the k-th. central moment can be obtained easily 
since terms with index t < k in (fTT|) will be asymptotically negligible. Indeed, we will use ®. 



Recurrence of M. k (s). We now focus on properties of AA k (s). By writing (flOl) in the form 

we see that -Mfc(s) satisfies the recurrence 

M *(*) = T~Pis) S (* > !)> (12) 

0<j<fc ^ 7 

with M (s) = 1/(1 - P(s)), where 

0<£<fc ^ 7 
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For example, Mx(s) = B 1 (s)/(1 - V{s)) 2 = (V(s) + pP'(s))/(l - V{s)) 2 . 

Note that each Bk(s) is analytic for 3?(s) > k and, in particular, for 9ft(s) > p. Moreover, the crucial 
property here is 

B x (p)=V{p)+i<P'{p) = 1-1 = 0, 

by our construction. Similarly, 

BaOO = W + 2^'(p) + pV(p) = p 2 P"(p) - 1 = a 2 /p, 

and B[{p) = a 2 /fi 2 . 

On the other hand, by (PT21) . we see that .Mfc(s) is analytic for 3?(s) > p and for 3?(s) = p except at 
s = p. Furthermore, because B\(p) = 0, it follows by induction from (PT21) . that at s = p, A^a,(s) has a 
pole of order at most [k/2\ + 1. 



Even moments. More precisely, for even k = 2£, we get by induction 

M k (s)~c k {s-p)- k / 2 -\ 

where 

k \..o/-\- k(k-l)_ 2 



Ck = \^ 2 j/J,B 2 (p)c k -2 = 2 ° Cfc ~ 2 ' 

with c = p, which is solved to be 

Cfc = plyl k\. 
We now apply Delange's Tauberian theorem and obtain 



E{Y N -p\ogv N ) k = — — h M 



n<N 
Ck 



N p (\ogN) k/2 



pT{k/2 + l)A(N) 

—f——a k (\ogN) k l 2 . (13) 
2 k / 2 (k/2)\ v 6 ; 



Odd moments. Let now /c = 2£ — 1, £ > 1. Since the coefficients fofc(n) are not necessarily nonnegative, 
we cannot directly apply Delange's Tauberian theorem. Instead, we consider the following two auxiliary 
Dirichlet series 

:=^ ra ~ S 5Z am W ((^-plogn) fe + (logn) fc/2 ) 2 , 

n>l m>0 

and, see © and ©, 

AjO) :=^n~ s ^a m (V) ((m-plogn) 2fc + (logn) fc ) 

n>l m>0 

= A<»(a) + (-l)MW(a). 
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The two Dirichlet series have only nonnegative coefficients, and we will show that Delange's Tauberian 
theorem can be applied to both series. The leading terms will cancel and we will have 

~JTpr\ E &*(n)0°g")* /a = ° ((^gN) k ) . (14) 

^ ' n<N 

From this, we use the monotonicity of (log n) k l 2 and elementary arguments to recover the desired estimate 

E(Y N -^\ogv N ) k = — L- h(n) = o ((logiV) fe / 2 ) . (15) 

^ ' n<N 

Proof of dH). Let 

D 3 (s) :=J2(^) k/2 bk(n)n- s . 

n>l 

Then -Di(s) = D 2 (s) + 2D 3 (s). By the discussions above, we can apply Delange's theorem to D 2 (s) and 
obtain 

^^^a^nl^-zilognf + flognf) ~C k (\ogN) k , (16) 

^ ' n<Nrn>0 

where C k = (2k)\a 2k / (2 k k\) + 1 (this value is however immaterial). 

We now show that the partial sum of the coefficients of Di(s) has asymptotically the same dominant 
term. We start from the representation (k = 2t — 1) 



D 3 (s) = (-lYn- 1 / 2 



POO 

Jo 



for 9fts > p, because (— l)fM^ f (s) = Yl n >i b k (n)(\ogn) e n s and 

poo poo 

(-l)W M[ l \s + t)r^ 2 o\t = y^b k (n)(\ognYn- s e -*iog« r i/a dt 
Jo ' n > x Jo 

= r(i)^6 fe H(logn) fe / 2 -- s 



n 

n>l 



We now consider the local behavior of D 3 (s) when s ~ p. First, .Mfc(s) has a pole at s = p with leading 
term c' k (s — p)~~( k+l ^ 2 , for some c' k . Thus A^^(s) has a pole with local behavior c'{,(s — p)~ fc_1 . It follows 
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that for small \w\ and 3?(w) > 0, 

POO 

D 3 (p + w) = {-If-R- 1 ' 2 / Mf (p + w + t)r^ 2 dt 

Jo 

= o(J™\w + t\- k -H-V 2 dtj 



r\w\ poo 

O [ / \ w \- k - 1 r 1 ' 2 dt + / r fc - 3/2 dt 

'0 J\w\ 



0(\w\- k - 1 ' 2 ) 
o(\w\- k ~ 1 ). 



Since Z?i(s) = D 2 (s) + D 3 (s) has all coefficients nonnegative, we can now apply Delange's theorem to 
D 1 (s) and conclude that 

— X) E (( m - ^ lo S n ) fe + ( lo g^) fe/2 ) 2 ~ Cfc(logiV) fc . 

^ ' n<Nm>0 

This, together with (TT6b . proves (PT4l) . 
Proof of (HU). Let 

B k (x) :=^bk(n)(\ogn) k / 2 , 

n<x 

and for a; > 2, cp(x) := (loga:)~ fe / 2 . Then 

/N t-N 
B k (x)ip'(x) dx = ^2 h(n)(\ogn) k/2 ip'(x)dx 
2<n<N ^ n 

= Y, ^H(logn) fc / 2 ((logiV)^/ 2 -(logn)- fc / 2 ) 

2<n<N 

= B k {N){\ogN)- k ' 2 - h(n). 

2<n<N 

Thus, by f[T4]>. 

-AT 



Kn<JV 



^ b k {n) = b k {l) + B k {N){\ogN)- k ' 2 - B k (x)<p'(x)dx 

= 0(l) + o(^(logAT) fc / 2 ) + - / B k (x)x-\\ogx)- k/2 - l dx 

= 0(1) + o (N p {\ogN) k ' 2 ) + o a^-^loga;)*/ 2 - 1 dz 
= o(i\^ (log iV) fc/2 ) , 



as required. 



From Yn — \i log v N to Y N — fi log N. The two estimates (1131) and (1151) imply 

/iWbg^y J jtf^rv iffciseven (17) 

v ^vioiiv ; ^ / if fc i SOdd , 



which in turn implies, by the method of moments, that 

Y N - ji log v N d 



^(0,i; 



av^ogiV 

Our final task is to prove the same asymptotics © from the two estimates (fT3l and (fT5l) . To that 
purpose, define S^x) = if x < 2 and 

S fc (x) := ^6 fc (n) (x>2). 

n<a; 

We use © and the cruder estimates (by ([[3]), (Q~5]) and ©) 

= O (A(x)(logx)' /2 ) = O (x p (\ogx) e/2 ) , 
for £ = 0, . . . , k — 1. With this, we have 



Now for each £ = 0, . . . , k — 1, 



jT (log ^ d$(z) = O uf " (log ^ aT 1 ^) dx 



Splitting the integral at x = N/2, and making the change of variables x 1— > iV/a; for the first half, we see 
that 

»iV / j\r\ k-e-i 



J (log— J x p - l {\ogxf' 2 dx 

/ rN/2 

= 0\ N p / x-P-^Xogxf- 1 - 1 



( N\ e/2 r N \ 

^log — J dx + J x^^ogx)^ 2 dx j 



= 0(N P (log N)^ 2 ) 
= o(N p [log N) k / 2 ) , 
for < £ < k — 1. This proves that 



Ei 



(Y N -^\ogN) k = h(n) + o ((logiV) fc / 2 ) 



n<N 

and thus the estimates in © hold by CE3) and (fT51) . 
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An alternative argument. To bridge (flTT) and we can also argue as follows. Consider the sum 

1 / N\ k 

E(\ogN-\ogu N ) k = —— £a(n) log- , 

^ ' n<N ^ ' 

which is 0(1) by a similar summation by parts argument as used above. Then for even k 

IJlogiV — logi/jv-IU = O(l). 

By Holder's inequality, this holds true also for every k > 0. Consequently, using again Holder's inequality, 
we deduce (Q~|). 

3 Conclusions and additional remarks 

While a direct application of Tauberian theorems leads to results of the form 

E(Y N ) -/xlogiV, 
we indeed prove, still relying on Tauberian arguments, that 

E(Xn) = filogN + o(y/fo£N), 

(a special case of Proposition d]). This shows the power of our approach. However, the estimates ([I]) 
we derived are not strong enough so as to prove more effective bounds such as the convergence rate to 
normality (or the Berry-Esseen bound). 

Another corollary to our moment convergence result is the following asymptotic approximations to all 
absolute central moments 



E\Y N -/ilogNf ~ 2 /3 / 2 7r~ 1 / 2 r ( 



(3 + 1 



(logiV) 



13/2 



for all (3 > 0, which seem difficult to get directly from Dirichlet series. 

On the other hand, when z E (— \ogV(n), oo), one can apply directly Delange's Tauberian theorem to 
the generating function 

1 

1 - e z P(sY 

(instead of to the Dirichlet series of higher moments obtained above by Taylor expansions in z); this results 
in the asymptotic approximation 

E (e zYN ) pP '^ NM-p 

where p(z) solves the equation 1 = e z P(p(z)) with p(0) = p. From this approximation, one might expect 
asymptotic normality by straightforward argument. However, the asymptotic result so obtained holds only 
pointwise, and the uniformity in z is missing here. While the gap of uniformity may perhaps be filled 
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by applying suitable Tauberian theorems with remainders, the use of Delange's Tauberian theorems is 
computationally simpler and technically less involved. 

It is clear from our proof that Theorem \T\ actually holds for any Dirichlet series V(s) with nonnegative 
coefficients and satisfying the conditions of Theorem [TJ Thus the restriction of 8? to a subset of positive 
integers is not essential. For example, one can consider the ordered totient factorizations with V(s) = 
E n >3^( n )" S ' where 0(n) is Euler's totient function, namely, the number of positive integers < n and 
relatively prime to n. In this case, k — 1 and p m 2.26386 since 

5>(n)- = coo n (i-p- s + ( P -i)- s ); 

n>l p: prime 

see for a detailed studied of this Dirichlet series. 

How to compute p to high degree of precision? In general, the problem is not easy, for example, if 
^( s ) = XL>2 r n/3 Q°S n ) c l ~ s f° r A c > 0' tne case °f totient factorization is similar. The easy cases 
are when V(s) can be expressed in terms of ^-functions such as V(s) = ((s) — 1 (all integers > 1) or 
V{s) = ((s)/((2s) — 1 (square-free integers > 1). Take now V{s) = J2 p . • e p~ s . The zero ofV(s) = 1 
can also be easily computed by using the relation (see tfTTI ) 

p : prime fc>l 

where p(k) is Mobius function. This readily gives 

p w 1.39943 33287 26330 31820 28072 

and 

p w 0.57764 86251 95138 05440 61351 
a 2 w 0.48439 65045 13598 28128 07456 .... 

We indicated a few directions to which our approach can be extended in Introduction. But can a similar 
idea be modified so as to deal with arithmetic functions with mean other than log TV? 
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